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Abstract 



This work deals with the computation of electron pair correction to 
small angle Bhabha scattering, in order to contribute to the improve- 
ment of luminometry precision at lep/slc below 0.1% theoretical 
accuracy. 

The exact QED four-fermion matrix element for e + e~ — > e + e~e + e~ , 
• including all diagrams and mass terms, is computed and different 

Feynman graph topologies are studied to quantify the error of ap- 
proximate calculations present in the literature. 

Several numerical results, obtained by a Monte Carlo program with 
full matrix element, initial-state radiation via collinear structure func- 
tions, and realistic event selections, are shown and critically compared 
with the existing ones. 

The present calculation, together with recent progress in the sector 
of 0(a 2 ) purely photonic corrections, contributes to achieve a total 
theoretical error in luminometry at the 0.05% level, close to the current 
experimental precision and important in view of the final analysis of 
the electroweak precision data. 

Keywords: electron-positron collision, small angle Bhabha scat- 
tering, theoretical error, light pairs, Monte Carlo 
PACS: 02.70.Lq,12.15.Lk,13.40.Ks,13.85.Hd 
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1 Introduction 



The high-precision determination of the machine luminosity at lep / slc is 
an essential ingredient of the success of precision tests of the electroweak 
interactions on top of the Z resonance 

As well known, the Bhabha scattering process at small angle (of the or- 
der of a few degrees) is the reference reaction used for luminosity monitoring 
at lep/slc, owing to its large cross section (dominated by t-channel pho- 
ton exchange) and its substantial independence of purely electroweak effects. 
Experimental efforts in the development of efficient, dedicated luminometry 
detectors, as well as precision calculations of the small-angle Bhabha (here- 
after sabh) scattering cross section both contribute to achieve a measurement 
of the lt Z factories" luminosity with a total relative error at the 0.1% level 
H], 0, [| . On the experimental side, the present total uncertainty is smaller 
than 0.1% @, close to the 0.05 level ||. As far as the theory contribution 
to the luminosity measurement is concerned, the estimate of the theoretical 
errors, used by the lep collaborations, is summarized in table [l] || for centre 
of mass energies around and above the Z resonance. 

Table 1: Theoretical error in sabh scattering according to ref. H at typical 
LEPl and lep 2 energies. 



Type of correction/error 


LEPl (%) 


LEP2 (%) 


missing photonic 0(a 2 L) 


0.100 


0.200 


missing photonic 0(a 3 L 3 ) 


0.015 


0.030 


vacuum polarization 


0.040 


0.100 


light pairs 


0.030 


0.050 


Z-exchange 


0.015 


0.000 


total 


0.110 


0.250 



Some comments on table [l] are in order. The components of the theo- 
retical error refer to the sabh scattering cross section, for any typical event 
selection of LEP experiments, as computed by the program BHLUMI v4.03 
H. The largely dominating source of theoretical error is due to the missing 
part of 0(a 2 L) subleading photonic corrections, where L = ln(— t/m?) is 
the collinear logarithm in t-channel scattering. Also the contribution of the 
missing part of the leading 0(a 3 L 3 ) corrections is of photonic nature. The 
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vacuum polarization entry is the effect of the uncertainty in the hadronic 
contribution to the running of «qed, when considering the parameterization 
and relative error estimate of ref. 0. The next contribution is the uncer- 
tainty introduced by the corrections due to the production of light pairs, 
chiefly e + e~ ones. The last entry refers to the uncertainty associated to the 
treatment of the 7-Z interference. More details about the strategy adopted 
in order to estimate the various sources of theoretical error can be found in 
ref. 1. 

After the analysis of ref. ||, important theoretical developments took 
place. Additional work in the sector of two-loop photonic corrections J?]. U 
led to the conclusion that the perturbative contribution due to the uncon- 
trolled part of 0(a 2 L) corrections does not exceed the 0.03% level. This con- 
clusion has been very recently reinforced by the detailed analysis of ref. Q. 
Furthermore, new determinations [0] of «qed lower the error on hadronic 
contribution to vacuum polarization in s-channel processes at y/s = Mz- 
This might affect sabh scattering too, although no dedicated analysis for 
the low- angle regime exists yet. 

As a consequence of this progress, it is relevant to reduce the uncertainty 
associated to the light pair contribution. At present, the calculations avail- 
able in the literature and used to estimate the light pair uncertainty as given 
in table [I] concern a Monte Carlo (hereafter mc) computation based on an 
approximate t-channel matrix element [[Hj and to an analytical approach 
with fixed event selections |TIJ|. Previous leading logarithmic evaluations of 



the dominant light pair contribution to SABH can be found in ref. |L2 . 

In order to improve the existing situation and contribute to the lower- 
ing of the light pair error, in the present paper a MC calculation is drawn 
with the exact e + e~ — > e + e~e + e~ matrix element and taking into account 
realistic event selections. The analysis is chiefly presented at lepI/slc en- 
ergies (v^i = 92 GeV), but numerical results are shown at lep2 energies 
(y/s = 176 GeV) too. The impact of the present calculation in the reduction 
of the theoretical error for lep / slc luminosity measurement is also discussed. 

The outline of the paper is as follows. The details concerning the treat- 
ment of phase space are described in section [|, while in section ^ the selection 
criteria considered in the present study are reviewed. Section f| is devoted to 
describe the calculation. The last sections contain a discussion of numerical 
results, including comparison with existing analytical calculations (section 
H) and approximate MC results (section ||), as well as study of the effect 
of initial-state radiation (hereafter isr) (section [?]). The conclusions and 
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possible developments are given in section 



2 Phase Space Generation 

Pair corrections to Bhabha scattering lead to a four-body kinematics and to 
an 8-dimensional phase space. Expansion in few body processes can greatly 
simplify the phase space parameterization. The choice between different 
equivalent expansions should be suggested by the relevant dynamics. 

At high energies and small momentum transfer, which is of interest for 
luminosity measurements at lep, the leading contribution to cross section 
is given by bremsstrahlung Feynman graphs sketched in figure [I], i.e. by 
t-channel photon exchange dynamics. Indeed, since the relevant dynamics 
is dominated by pure QED processes, weak effects, such as Z-exchange or 
multiperipheral graphs mediated by at least one Z boson, will be neglected 
in the following. Further the present study deals mainly with electron pair 



production, because it is known |3|, [TT| that heavier particles give a much 
smaller contribution. 



P- qi 




p+ qi 

Figure 1: One of the sixteen bremsstrahlung graphs representing the leading 
t-channel dynamics. 

Let us define the photon momentum = g 3 + g 4 and its energy by to = ko. 
The core of bremsstrahlung contribution is given by the soft pair approxi- 
mation, i.e. the limit \t\ 3> u, \k\. In this regime the emitted pair is almost 
collinear to the photon k. Thus the phase space configurations in which q 3 
and ^4 are back-to-back are highly suppressed by t-channel dynamics. 

However, the selection criteria for kinematic events, used by the lep 
collaborations and reviewed in section |3|, scan also the hard region. When 
bremsstrahlung processes get smaller, the next to leading Feynman graph 



3 



topology is represented by multiperipheral dynamics shown in figure |2|. No- 
tice that this contribution is relevant also for 77 physics, being described in 
its bulk by the Weizsacker- Williams approximation ||13| for which the internal 
photons become quasi-real. 

P- qi 




p+ q-i 

Figure 2: One of the eight Feynman diagrams for multiperipheral dynamics. 

Bremsstrahlung and multiperipheral graphs do not complete all the Feyn- 
man graph topologies. Other two classes of diagrams can be drawn, namely 
the annihilation and conversion ones, which are shown in figure |3[ Their con- 
tribution is less important at high energies and small momentum transfer. 
Thus in this paper phase space parameterization and importance sampling 
does not deal with these configurations. 

P+ qi P- qi 




P- q-2 P+ q-2 

Figure 3: Two of the twelve Feynman diagrams representing conversion and 
annihilation dynamics, respectively. 

The two following subsections show how the kinematics is treated accord- 
ing to the previous considerations about the dynamics. 

2.1 Phase Space for Bremsstrahlung Events 
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2.1.1 Phase Space Parameterization 

Let us consider the four-body phase space for two incoming particles of mo- 
mentum p_ and p + , and four outgoing particles of momentum q\ . . . 54. Let 
us define P = p_ + p + . 

d s R 4 (P; Ql , q 2 , q 3 , q 4 ) = {[ ^ 4 \P - £ ft ) (1) 

i=l 1 i=l 

The four-body kinematics can be split into the product of simpler pro- 
cesses, t-channel dynamics suggests to think of the event as formed by a 
fermionic current interacting with a particle which radiates a pair. Thus let 
us deal with the event as a three body plus a two body kinematics. The 
following decomposition can be introduced 



d 8 R 4 (P; q u q 2 , q 3 , q±) = d 5 R 3 (P; q u q 2 , qu)d 2 R 2 (q u ; 53, Qi)ds 34 (2) 

where s 34 is the squared mass of g 34 = q 3 + g 4 . 

The three body kinematics can be expressed in terms of invariants. Let 
us define s^... = + qj + . . .) 2 and t±j,„ = {p± — qj — . . .) 2 . A straightforward 
calculation leads to 



d 5 R 3 (P; qi, q 2 , qu) = 00 . 1/2 , = 2 , ^-r- ds 12 dt^ 34 ds lu dt +2 d(p (3) 

where dy? integration is over the beam direction, A is the Kallen function and 
A 4 is the fourth order symmetric Gram determinant. 
The two body kinematics gets the simpler form 

, 2p/ x X 1/2 (s 3A ,m 2 ,m 2 ) 2 

d R 2 {q 3 4, 5s, 54) = 5 d (4) 

oS 34 

where f2 3 is the solid angle of q 3 in the frame in which <f 34 = 0. 
2.1.2 Physical Region 

The physical region depends on the choice of the order of integration for the 
phase space variables. Let us consider the following choice (read up-down 
then left-right) 
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(5) 



These limits are exact, thus the phase space generation has efficiency equal 
to one if no selection criterion is set on. 



2.1.3 Importance Sampling 

The soft pair limit provides a relatively simple analytical approximation 
for the t-channel contribution to the cross section. This expression can be 
reached either from direct calculation or from analytic results available in 



the literature |L4] . The resulting integral gives a guideline to sample the full 
cross section formula. 

According to the choice for the phase space variables ordering, the weights 
w follow 
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w{s 34 ) 
w(s 12 ) 



1 

S34 

f 1/(S + S34 - S12) OC (I/E34) ,Sl2>si 2 

1 1 , S12 < s' 12 



^(*-34) = -5-^— (6) 



m 2 — t 



34 

1 



S134 - ra 2 



1- , t+2 < *+2 

w(t +2 ) = <| l/(*2 2 ) ,t' +2 <t +2 <t% 
l/(-t +2 ) ,t +2 >t% 

where w(s 34 ) and w(si 2 ) deal with the infrared pole, w(t_ 3 4) and w(si 3 4) 
with the collinear pole, and w(t +2 ) with the coulomb pole. The boundaries 
s' 12 , t' +2 and t" +2 should be set according to the selection criteria. 



2.2 Phase Space for Multiperipheral Events 

Multiperipheral events require true four-body kinematics, since no natural 
expansion of the phase space can fit the dynamics, because of the three prop- 
agators between the ingoing particles. Thus a description based on energies 
and angles can be preferred to an invariant picture. 



2.2.1 Phase Space Parameterization 

Let us integrate out g 4 and put E; L = q i 

<eUP; «.,*,*,«) = HE - ± ^g^g^L (T) 

where Ei, E 2 and E 3 are given by mass shell relations and E4 by energy- 
momentum conservation. 

Let us consider the decay k — > ^3,^4 where k = q 3 + g 4 , u = k , and 
s 34 = k 2 . The three momentum |<f 3 | and the cosine of the polar angle 9 3 , 
with respect to a fixed direction u z , are needed as phase space variables. Let 
us write k as k = k z u z + k x il x the third axis being defined as u y = u z x u x . 
Energy conservation implies that 
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namely 



uj 2 + E% - 2uE 3 = m 2 + \q 3 \ 2 + \k\ 2 - 2k z \q 3 \ cos0 3 + 

-2k x \q 3 \sm9 3 smip 3 (8) 



siny9 3 = —————— (2a; £"3 - s 34 - 2k z \q 3 \ cos0 3 ) (9) 

2k x \q 3 \ sin 3 



The Jacobian is easily obtained as 



j-i = k x \g 3 \ sin 3 cos y? 3 
£" 4 



This leads to the following phase space parameterization 



d s R A (P]q l ,q 2 ,q 3 ,qA) = — . x 

lefc^i^-Eal&l sm 3 cos y? 3 

<i|gi|(i 2 f2i(i|g2|rf 2 ^2rf|93Mcos03 (11) 
where k x — (P — q\ — q?) • u x , Ei = ^m 2 + |$|, and <i 2 f2j = d cos 9idifii. 

2.2.2 Physical Region 

The physical region is given by 

|2/c x |93| sin^sl > ^o^Ea — uj 2 + s 34 — 2& 2 |(f 3 | cos0 3 | (12) 
which means (sin 3 > 0) 



2{}z x \q 3 \ sin 3 - u E 3 + k z \ q 3 1 cos 3 ) > -u; 2 + s 34 
2(/c x |g 3 | sin 3 + uj E 3 - k z \ q 3 | cos 3 ) > uj 2 - s 34 



(13) 



By taking into account that \q 3 \ < \J\q 3 \ 2 + m? < \q 3 \ + m it follows that the 
range of interest is certainly contained into the domain 

2(fc x sin0 3 - u + k z cos9 3 )\q 3 \ > -u? + s 34 . . 

2{k x sn\9 3 + uj — k z cos9 3 )\q 3 \ > uj 2 — s 34 — 2ujm 



which leads to 
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2(k x sin # 3 + lo — k z cos # 3 ) ~~ ~~ 2{oj — k x sin # 3 — k z cos #3) 
if a; — sin # 3 — A; 2 cos #3 > 

The maximum and minimum of |g*3| are achieved for <f 3 collinear to k in the 
forward and backward direction, respectively. 

2.2.3 Importance Sampling 



The Weizsacker- Williams equivalent photon approximation [13] provides a 
useful guideline to deal with the contribution of multiperipheral diagrams. 
The resulting integral shows how to sample the full cross section formula by 
introducing the weights 

1 1 

w(cos#i) = — — — , w(cos6 2 ) = — — r- (16) 

-i_i(cos6>i) — i +2 (cos0 2 ) 

where w (cos ^1) and w(cos 82) mimic the most singular behaviour of the Weiz- 
sacker- Williams spectrum. t_i and t +2 are equal to (p_ — qi) 2 and (p + — g 2 ) 2 , 
respectively, as previously defined. Further a flat importance sampling can 
be made on cos #3 to match the selection criteria better. 



3 Selection Criteria 

Two types of selection criteria are considered: BARE1 and CAL02. They were 
defined by the "Event Generators for Bhabha Scattering" working group 
during the cern Workshop "Physics at lep2" (1994/1995), see ref. [|J for 
more details. 

These algorithms are tailored to Bhabha scattering with photon emission. 
Thus they must be modified to fit with pair emission, chiefly to deal with 
identical particles and with fermion clusters. 

3.1 A Simple Setup: BARE1 

This non-calorimetric criterion selects events with one or more particles per 
calo. If more than one particle hits a calo it chooses the most energetic 
hit, the reference particle. Thus for each kinematic event two particles are 
labelled as the emitted pair and the other two as the beam particles. 
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Table 2: Angular acceptances in degrees for BARE1 and CAL02 algorithms. 





BARE1 


CAL02 


Left Calo 


Right Calo 


Left Calo 


Right Calo 


WW 


2.70 -> 7.00 


2.70 -> 7.00 


2.97^ 6.73 


2.97^ 6.73 


NN 


3.30 -> 6.30 


3.30 -> 6.30 


3.49 -> 6.11 


3.49 -> 6.11 


NW 


3.30 -> 6.30 


2.70 -> 7.00 


3.49 -> 6.11 


2.97 -> 6.73 



Then the beam particles are selected by their energies and angles. The 
angular cut can be made with the same acceptance for the two caloes (wide- 
wide and narrow-narrow), or with different acceptances (narrow- wide) . 

3.3° < 6 < 6.3° — narrow acceptance , , 

2.7° < 9 < 7.0° - wide acceptance ^ ' 

Angular acceptances are summarized in table 0. Further the energy cut is 
imposed by considering the adimensional variable z as usually done in the 
literature JI£| 

z = l- — 2 — < z max (18) 

beam 

It is worth noticing that small values of z means soft pair emission, while 
large values allows for hard pair radiation. 



3.2 A Calorimetric Setup: CALQ2 

This selection criterion looks for any particle in two caloes, the former is 
along the ingoing electron direction, the latter is along the ingoing positron 
direction. Then it scans the detected particles in each calo to pick up the 
most energetic one, the reference particle. 

It defines a cluster for each reference particle. The shape of the cluster is 
a square in the (6, ip) plane of the corresponding calo, its size is of 3A(9/16 x 
37r/16 radiants, where A8 is the calo polar width. The cluster center is pinned 
to the reference particle. The cluster energy is the sum of the energies of the 
particles in the cluster itself. 

The last step of the algorithm is to reject the events according to an 
angular and an energy cut, where angles are referred to the cluster center 
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Backward hemisphere 



WW: Of G (C in ,CaJ, NN: Of £ (C..C 
NW: 6? e # 6 (C-«L) 



2tt 



Figure 4: 
ref. [HI. 



Geometry and acceptances of CAL02 calorimetric setup. From 



and energies to the cluster energy. Each polar width must be reduced to an 
effective acceptance, as shown in figure |], to prevent the cluster square to 
exceed the calo size, i.e. should the reference particle be out of the effective 
acceptance but within the calo total width, yet the event is rejected. 

The angular cut can be made with the same polar width for the two 
caloes, defining a wide-wide and a narrow-narrow acceptance, or with differ- 
ent widths for the two caloes, defining a narrow-wide acceptance 



3.3° < 6 < 6.3° - 
2.7° < 6 < 7.0° - 



width for narrow acceptance 
width for wide acceptance 



(19) 



Angular acceptances are summarized in table [| Further the energy cut is 
imposed by considering the adimensional variable z as given by eq. (ffH). 



4 Dynamics Calculation 



Generators including the full set, or part, of the QED diagrams for e + e — > 
(Z = e, /i, r) was already described in ref. [0 and they are used 



for analysis and simulation of 77 collision processes at lep and other e + e 
colliders. 
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In the present paper the MC program PAIRS, written in FORTRAN, was built 
to compute the e + e~ — > e + e~e + e~ process. It implements the importance 
samplings and the selection criteria sketched previously, and computes the 
exact matrix element, including all mass terms, by using the ALPHA algorithm 
and the resulting code |15|], that is conceived for the automatic computation 
of tree-level multi-particle production amplitudes without any need of Feyn- 
man graphs expansion. 

Without entering the details of the algorithm ALPHA it is worth noticing, 
for the aim of the present study, that the predictions of this automatic algo- 
rithm have been already compared with diagrammatic results for processes 
with four fermions in the final state [17] showing an excellent agreement, and 
also successfully used to obtain original results for other reactions 



The main features of the program PAIRS can be summarized as follows. 
The code computes the phase space integral by means of an importance 
sampling both for bremsstrahlung and multiperipheral graphs. Since the 
contribution of the other topologies is small, as discussed in section |2|, there 
is no need of a specific strategy to reduce the associated variance. 

The integration is performed in two steps. In the former the matrix 
element is sampled by using the bremsstrahlung weights, see section |2. 1.3 . 
Then a rejection algorithm selects the events belonging to the phase space 
region Qb which is preferred by the bremsstrahlung dynamics. In the latter 
the same calculation scheme is performed in terms of the multiperipheral 
weights, see section |2.2.3| , and of a rejection algorithm selecting the remaining 
phase space region Om, which is preferred by the multiperipheral dynamics. 
The treatment of the interference between the two dynamics is addressed 
below. 

After phase space generation the program deals with identical particles 
and other symmetries of the amplitude. Feynman expansion in pure QED 
shows sixteen bremsstrahlung graphs, eight with initial-state emission and 
eight with final-state emission of a fermion pair. These two families can be 
sampled at the same time leading only to eight different weights, which can 
be obtained, fixed one of them, by repeated application of identical fermion 
exchange symmetry (ID symmetry) or of CP symmetry. Let us note that 
under these symmetries the squared amplitude does not change. There are 
also eight multiperipheral graphs, four with an amplitude value and four, 
obtained twisting the kernel fermion lines (TW symmetry), with a different 
value. The two multiperipheral squared amplitudes are left invariant by the 
ID symmetry. 
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The symmetries of the integrand allow us to cast the cross section for- 
mula in a simpler form, because the phase space jacobians and the selection 
criterion characteristic functions share these symmetries too. 

An example can clarify the importance of the integrand symmetries. Let 
us consider a function f(x,y), symmetric under the change x — * y, which 
is integrated over the unit square in the (x, y) plane. In order to draw the 
calculation some weights are introduced, for defmiteness say iO\ = l/(x + a) 
(a > 0) is a suited weight because of some pathology in x = —a, but, since 
f(x, y) is symmetric under coordinate exchange, it must also exhibit the same 
pathology in y — —a. Thus the weight u>2 = l/(y + a) must be introduced 
too. The integration by using importance sampling follows. 

f 1 dx f 4y f( X ,y) = f 1 dx f 4y ( + - 

JO JO Jo Jo \ iO-\ + LUo 



o x + a Jo x + y + 2a 

+ [ 1 dx[ 1 ^ ( X + a) ( y + a) f(x 1 y) = 
Jo Jo y + a x + y + 2a 

1 /-ln(a+l) ( x + a \( y + a ) 

dx / dfi - " - f& y) 20 

Jino x + y + 2a 

In the last step coordinate exchange symmetry both of f(x, y) and integration 
region was used, and the integration measure changed in order to include the 
weight by putting dfi = dy/ (y + a). 

Hence the illustrated two channel sampling can be collapsed into a single 
integration, because the integrand and the integration region share the same 
symmetries. The same strategy can also be adopted for the more complex 
integral involved in the present calculation. 

Let us name by x the phase space coordinates, by d 8 R 4: (x) the phase space 
volume element defined by eq. (|]), by iV the cross section normalization, by 
x{x) the characteristic function of the selection criterion, by pf(x) the eight 
sets of weights for the bremsstrahlung amplitude and by pf 1 (x) and pf^xw) 
the four weights for each of the multiperipheral amplitudes. These quantities 
share the symmetry properties 

d 8 R 4 ( Xi ) = d 8 R 4 (x) \/ie{ID,CP,TW} 

X (xi) = X {x) \/iE{ID,CP,TW} 

\M{ Xl )\ 2 = \M{x)\ 2 Vie {ID, CP} 

pf' M ( Xj ) = pf M {xi) Vi,je{ID,CP,TW} 
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where phase space point obtained from x applying a suited symmetry 

among ID, CP, and TW. 

In terms of these quantities the cross section integral becomes 



a — sn [ dH ^ x)p (x ^vi^)i 2 + 

m j m^M*) E (|M(x)|2 + Mxtw)12) m 

l^ie{ID,TW}Pi \ x ) spin 

This technique is very useful because it permits to treat twenty-four channels 
as they were three. 

The above division of phase space leaves some problem in the brems- 
strahlung region because of the interferences between t-channel and multi- 
peripheral dynamics. Thus it is useful to add an extra flat channel to the 
first integral of eq. ([22|). The flat channel deals also with a multiperipheral 
contribution so it must be sampled further by the p^(x) weights. Hence the 
cross section integral is 



a 



sn[ El^)l 2 + 

8Nrj ( ^(*)r"(*)x(*) ^ E | M(X) | 2 + 

m j ^Kfe^ E (|M(x)|2 + mxTw)?) m 



spin 



where i] is a suitable weight for the flat channel. 

The matrix element is computed by functionally integrating the qed tree 
level effective lagrangian with the iterative algorithm ALPHA, which gives 
directly the total amplitude M(x). 

After the cross section is computed for the real contribution, also the 
virtual correction corresponding to pair emission is added. The resulting 
pair contribution is then normalized against the tree level Bhabha scattering 
cross section. The final result can be corrected for initial-state radiation via 



collinear structure functions too 20 . 
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5 Comparison between MC Results and An- 
alytical Calculations 

Before the numerical results being worked out, the MC program is now tested 
against analytical results already present in the literature Such calcu- 



lations are available for fixed selection criteria, among which one resembling 
BARE1 algorithm. 

A direct comparison, which could be affected by some bias because of 
analytical approximations and possibly small differences in the selection cri- 
terion, is shown in figure [5] and table |3] for the absolute value of the pair 
correction. Solid line represents analytical computation of ref. |lTJ, while 
markers show the results of the present MC calculation. The error bars are 
within the markers. Entry values are in pb, sum up both real and virtual 
part, and are computed for the BARE1 setup with symmetric angular accep- 
tance 1.375° < 6 < 3.323° at ^/s = 92.3 GeV. 



Q_ 

b 
< 



20 


-20 
-40 
-60 
-80 
100 
120 
140 
160 
180 



0.2 



0.4 



0.6 0. 



Figure 5: Comparison between MC integration of the exact matrix element 
(markers) and analytical computation of ref. |Tl[ (solid line) as a function of 
the energy cut z, as given by eq . (p~8|) . Entry values sum up real plus virtual 
cross sections, and are computed for BARE1 setup with symmetric acceptance 
1.375° < 6 < 3.323° at ^/i = 92.3 GeV. 

Table |3| shows that the relative difference between numerical and analytical 
results is within 15% for the experimentally relevant values of z 
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Table 3: Comparison between present numerical results and analytical ones 
of ref. [TTJ with the same settings of figure [S]. Entry values are in pb. In this 
setup Bhabha tree level cross section is a = 175426(14) pb. The error on cq 
is due to numerical integration. 





z = 0.1 


z = 0.3 


z = 0.5 


z = 0.7 


z = 0.9 


Present MC 


-162 ±4 


-95 ±1 


-60 ±2 


-40 ± 3 


-13 ±3 


Analytical 


-142 


-90 


-58 


-33 


+7 



(0.3 < z < 0.7), thus providing a rather satisfactory test of the MC pro- 



gram. Further checks were performed with the results of ref. |TT]| for CAL02 
event selection leading to a similar behaviour. 

6 Comparison between Exact Matrix Ele- 
ment and ^-Channel Approximation 

Numerical computations already present in the literature are limited to brems- 



strahlung graphs without fermion exchange and up-down interference [10 
Let us define such an approach as t-channel approximation. This approxima- 
tion, together with the analytical work of ref. |T1[, is an essential ingredient 



to estimate the theoretical error associated with light pair correction to sabh 
scattering. 



As far as the MC calculation of ref. |10[ is concerned, the uncertainty 
is the sum of a physical error, due to the incomplete matrix element cal- 
culation, and a technical error, due to algorithm stability and finite CPU 
time. In ref. |10| the physical error associated to the matrix element for 



real pair production is estimated to be 30%. Once the virtual correction is 
added, the cancellation against the real part can increase the magnitude of 
the error. In ref. [|1(]] summing up all the sources of physical and technical 



uncertainty leads to a conservative estimate of the light pair correction of 
— 1.3 x 10~ 4 ±2 x 10 -4 for narrow-wide calo2 acceptance at z — 0.5 at 
lepI/slc energies. Thus the comparison of t-channel approximation with 
the full matrix element computation can size the magnitude of the physical 
error, and therefore lowers the theoretical error. 

In order to remove any spurious effect in the comparison between such cal- 
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Figure 6: In the upper row: exact matrix element (markers) vs t-channel ap- 
proximation (solid line) as a function of the energy cut z, as given by eq.(|l8D. 
Entry values sum up real plus virtual cross sections, and are computed by 
varying the angular acceptances of CAL02 setup at lepI/slc energies. In 
the lower row: relative difference between the exact matrix element and t- 
channel approximation still as a function of the energy cut. Entry values 
refer only to the real pair production cross section. 

culations and the present one, the amplitude has been directly calculated in 
t-channel approximation by using the ALPHA algorithm^ and comparing that 
amplitude with the full amplitude within the same selection criteria. Further 
the same virtual correction |Tj| for both the cross sections is implemented. 
The results are shown in figure |6| for CAL02 setup and different angular ac- 

1 Since the ALPHA approach does not make use of Feynman graph expansion, this 
is not entirely straightforward. The goal is achieved as follows: the input lagrangian is 
modified by introducing two distinct photons 71 and 72. The electron couples to both 
photons, whereas the muon and the tau lepton couples to 71 and 72, respectively. Within 
this modified lagrangian m T = m M = m e is assigned and the process e + /i~ — > e + /i~T + r~ 
is studied. It is immediately seen that this reproduces the t-channel approximation. 
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ceptances at lepI/slc energies. In the first row the solid line represents the 
t-channel approximation, while the markers the results of the present full cal- 
culation. Entries are normalized to Bhabha tree level cross section (<To) and 
sum up both real and virtual part. In the second row the markers represent 
the relative difference between the real part of the t-channel approximation 
(erf 1 ) and the real part of the full calculation (cr^ t ). 

The comparison shows appreciable differences only for z > 0.7, without 
much sensitivity to angular acceptance variation. By increasing z greater 
than 0.7, the relative difference between the real part of the cross sections 
can grow up to 15%. 

Table 4: Comparison between exact and t-channel approximation numerical 
results for narrow-wide CAL02 setup, see figure |7|. Entry values are in pb 
except for the last column. In this setup Bhabha tree level cross section is 
o"o = 21939(1) pb. The error on <jq is due to numerical integration. 





Full Dynamics 


t- Channel 


Abs.Diff. 


Rel.Diff. (%) 


z = 0.1 


-25.36 ±0.01 


-25.48 ±0.02 


0.12 ±0.02 


0.48 ±0.08 


z = 0.3 


-12.85 ±0.05 


-13.34 ±0.02 


0.49 ±0.05 


3.67 ±0.38 


z = 0.5 


-7.14 ±0.05 


-8.43 ±0.02 


1.29 ±0.05 


15.30 ±0.63 


z = 0.7 


-4.98 ±0.12 


-7.06 ±0.06 


2.08 ±0.14 


29.46 ±2.23 


z = 0.9 


-1.75 ±0.17 


-6.90 ±0.07 


5.15 ±0.19 


74.64 ±3.51 



A more refined computation is shown in figure [7|, allowing a higher statis- 
tics, to get an improved accuracy in order to establish a better evaluation 
of the physical error. The result is shown in figure [7] and table |] for CAL02 
narrow-wide setup. The solid line represents t-channel approximation, while 
the markers the full calculation. The error bars are within the markers. En- 
try values are in pb, and sum up both real and virtual part. Table [| shows 
the numerical values for the two calculations and the absolute and relative 
difference. Entry values are in pb. 

The relative difference between the t-channel approximation and the full 
calculation listed in table [| shows that the physical error, introduced by 
using the t-channel approximation, is within 30% for the relevant values of 
z. In the soft region (small z values) the relative difference gets smaller 
since the bremsstrahlung diagrams are the largely dominating contribution 
at this regime. On the other hand, in the hard region (large z values) the 
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Figure 7: Results of a refined computation of the exact matrix element 
(markers) vs t-channel approximation (solid line) as a function of the en- 
ergy cut z, as given by eq.(|18). Entry values sum up real plus virtual cross 



sections, and are computed for narrow-wide CAL02 setup. 



multiperipheral diagrams become more and more relevant, worsening the 
agreement with the t-channel approximation. 

Some tests performed at lep2 energies show the physical error behaviour 
given in table |5|. The difference between the exact calculation and the t- 
channel approximation is within 10% for the experimentally relevant z values. 

The previous results refer to the CAL02 calorimetric setup explained in 
section || where the calo acceptances are set according to ref. Hl0| . Yet 
some tests were drawn also with different angular acceptances to match the 
recent angular standard of lep luminometers, as previously considered in 
refs. 0, |J. The new ranges are (1.50,3.20) degrees for the wide acceptance 
and (1.62,2.84) degrees for the narrow acceptance at lepI/slc energies. 
The results are shown in table |6| confirming the same behaviour previously 
obtained. 

It must be emphasized that now a MC program, that includes all the QED 
Feynman diagrams for e + e _ — > e + e~l + l~ (/ = e, /i,r), is available to calcu- 
late pair production corrections. Hence the physical error for electron real 
pair production can be reduced well below the 30% quoted in the literature. 
Further it is worth noticing that muon pair production was checked to be 
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Table 5: Comparison between exact and t-channel approximation numerical 
results for narrow- wide CAL02 setup at LEP2 energies. Entry values are in 
pb except for the last column. In this setup Bhabha tree level cross section 
is cr = 6087.0(3) pb. The error on cr is due to numerical integration. 





Full Dynamics 


i-Channel 


Abs.Diff. 


Rel.Diff. (%) 


z = 0.3 


-4.47 ±0.01 


-4.27 ±0.01 


-0.20 ±0.02 


-4.61 ±0.24 


2 = 0.5 


-2.91 ±0.01 


-2.68 ±0.01 


-0.23 ±0.02 


-8.70 ±0.62 


z = 0.7 


-2.43 ±0.01 


-2.28 ±0.01 


-0.15 ±0.02 


-6.62 ±0.64 



Table 6: Comparison between exact and t-channel approximation numerical 
results for narrow-wide CAL02 setup at LEPl/SLC energies with a different 
choice for the angular acceptances: wide (1.50°, 3.20°), narrow (1.62°, 2.84°). 
Entry values are in pb except for the last column. In this setup Bhabha tree 
level cross section is cr = 87203(4) pb. The error on o"o is due to numerical 
integration. 





Full Dynamics 


t- Channel 


Abs.Diff. 


Rel.Diff. (%) 


z = 0.3 


-41.01 ±0.30 


-43.62 ± 0.20 


2.61 ±0.36 


5.99 ±0.85 


z = 0.5 


-23.91 ±0.42 


-28.91 ±0.25 


5.31 ±0.49 


18.36 ± 1.84 


z = 0.7 


-16.99 ±0.69 


-24.99 ±0.25 


8.00 ±0.74 


32.02 ±3.27 



one order of magnitude smaller than the electron one. 



7 Initial-State Radiation 

Corrections by ISR can be evaluated by the aid of collinear qed structure 



functions in the non-singlet approximation |f2(J . Numerical computations for 
the previously discussed t-channel approximation are known in the litera- 
ture [10]. Moreover, corrections of the order of 0(a 3 L 3 ), as due to single 



bremsstrahlung to pair production, was analytically calculated in ref. [|IT . 

In figure [8] the comparison between the exact matrix element and t- 
channel approximation is made in the presence of ISR with the CAL02 setup in 
the narrow- wide acceptance at lepI/slc energies and with the same nota- 
tion previously used. An apex 7 means that the relative correction is dressed 
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Figure 8: First picture: exact matrix element (markers) vs t-channel ap- 
proximation (solid line) as a function of the energy cut z, as given by eq.fllBD, 
in the presence of ISR. Entry values sum up real plus virtual cross sections, 
and are computed for narrow-wide angular acceptances of CAL02 setup at 
lepI/slc energies. Second picture: relative difference between the exact 
matrix element and t-channel approximation till as a function of the energy 
cut in the presence of ISR. Entry values refer only to the real pair production 
cross section. Third picture: relative effect of ISR correction for exact matrix 
element and t-channel approximation. 

by ISR. 

As can been seen, the presence of ISR becomes more and more important 
in the soft region (small z values) up to modifying pair correction by a 25%, 
but it does not alter the analysis in tree level approximation discussed in the 
previous section. 

8 Conclusions 

The progress in the reduction of the sources of theoretical error to Bhabha 
cross section at small angle, driven by increased experimental accuracy, 
showed the need of reducing the uncertainty associated to the light pair 
contribution. 

Till now the theoretical error to sabh scattering due to pair production is 
evaluated by using MC results based on t-channel approximation and with the 
aid of analytical means. In order to improve the present situation, a new MC 
program for e + e~ — > e + e~e + e~ process was presented. It includes the exact 
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QED four-fermion matrix element, ISR in the collinear approximation, and 
realistic selection criteria. The results obtained by means of that program 
were compared with the best approximations present in the literature, either 
analytical or numerical. 

Two kinds of contribution form the theoretical error: the physical er- 
ror, due to matrix element approximation, and the technical error, due to 
numerical integration. The comparison between the exact matrix element 
calculation performed in the present paper and the t-channel approximation 
adopted in the literature showed that the physical error committed by using 
such approximation at lep 1 / slc energies is within 30% with the exception 
of the hard region. Indeed in such region it grows up to 75% of the correction, 
while in the soft and intermediate region it can be much lower. Thus the cal- 
culation drawn with the t-channel approximation is within the physical error 
quoted in the literature [ffD| in the experimentally relevant kinematic region. 
At lep 2 energies the physical error due to the approximation is smaller and 
it is within 10%. 

Concerning the present calculation, the physical error associated to elec- 
tron pairs amounts in neglecting Z-boson contributions to the tree level ma- 
trix element, in approximating ISR by the collinear structure functions, and 



in omitting pair correction exponentiation as for instance done in ref. (10 
Relatively to the whole correction, the physical error due to these sources can 
be estimated within 5%, while the technical error amounts only to a few per 
cent, leading to a conservative 10% of the correction for the total theoretical 
error. 

Further light pair production sums up the electron pair production with 
other contributions, e.g. muon pair production. Such effect amounts for 
muon pair production to a 10-20% of the electron contribution. However, 
since the present algorithm can compute also the effect of muon pairs, the 
total theoretical error can be estimated at 30% level of the whole correc- 
tion, as due to approximations in electron and muon pair corrections and to 
neglecting heavier pair contributions. 

A review for the contribution to the theoretical error at lepI/slc is 
shown in table [F[ In the first column are listed the values of table [I] || and 
in the second column the updated values to the latest results for photonic 
0(a 2 L) corrections ||, [7|, |J, while in the third column the values are fur- 
ther updated to the result of the present analysis, when taking into account 
electron and muon pair corrections and a total error estimate of 30%. 

Table [7] shows that the theoretical error for sabh scattering, because 
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Table 7: Theoretical error in SABH scattering at LEPl/ SLC. The first column 
(see also table |l|) refers to ref . , the second column takes into account the 
results of refs. [7|, the third column updates to the present analysis. 



Type of correction/error 


ref. 1 (%) 


ref. i (%) 


updated (%) 


missing photonic 0(a 2 L) 


0.100 


0.027 


0.027 


missing photonic 0(a 3 L 3 ) 


0.015 


0.015 


0.015 


vacuum polarization 


0.040 


0.040 


0.040 


light pairs 


0.030 


0.030 


0.010 


Z-exchange 


0.015 


0.015 


0.015 


total 


0.110 


0.061 


0.054 



of the reduction of the pair production error, is now close to 0.05% and 
so comparable with the experimental accuracy. The conclusions, drawn for 
lepI/slc, can also be extended to LEP2 energies and lead to the reduction 
of the theoretical error shown in table H. 



Table 8: Theoretical error in SABH scattering at LEP2. The first column 
(see also table [I]) refers to ref. , the second column takes into account the 
results of refs. [|], [7|, the third column updates to the present analysis. 



Type of correction/error 


ref. § (%) 


ref. | 


(%) 


updated (%) 


missing photonic 0(a 2 L) 
missing photonic 0(a 3 L 3 ) 
vacuum polarization 
light pairs 
Z-exchange 


0.200 
0.030 
0.100 
0.050 
0.000 


0.040 
0.030 
0.100 
0.050 
0.000 


0.040 
0.030 
0.100 
0.015 
0.000 


total 


0.250 


0.122 


0.113 



It is worth noticing that now the major contribution to the theoretical 
error is due to the hadronic uncertainty in the vacuum polarization correc- 
tion. Thus, if the new determinations f| of «qed(A^z) could be drawn to 
the 1 GeV mass scale too, the total error could decrease further. 

The present work was conceived to lower the theoretical error on SABH 
scattering due to light pair production. Yet, in order to manage the numerical 
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and theoretical questions, techniques of a wider range of applicability were 
implemented. Chiefly, to deal with multiperipheral dynamics, Weizsacker- 
Williams approximation played an important role as a guideline to develop 
a suited importance sampling and it led to a numerical recipe that can be 
further applied to perform phenomenological studies of processes involving 
particles lost in the beam pipe. 

Therefore, future developments could deal with the evaluation of pair 
correction to two fermion production large angle processes at present and fu- 
ture e + e~ colliders. Experimentally relevant processes are e + e~ — > e + e~, 
whose interest is for luminosity measurements at da$ne and nlc, and 
e + e~ — > hadrons in the context of precision studies of the electroweak in- 
teraction at lep and beyond. 

Other fields of interest concern background processes for WW physics 
and searches for new physics, as very recently addressed in ref. [211] for the 
single W production case, and two photon collision processes. 

All these studies will require the upgrade of the QED matrix elements 
considered in the present paper to the full electroweak ones. It is worth 
noticing that this generalization could be worked out quite directly by virtue 
of the automatic algorithm employed in the present approach. 



Acknowledgements 

M.Moretti was supported by a Marie Curie fellowship (tmr-erbfmbict 
971934). A.Pallavicini wishes to thank infn Sezione di Pavia for the support 
and for computer facilities. 



References 

[1] See for example G.Montagna, O.Nicrosini, F.Piccinini, Riv.Nuovo Ci- 
mento 21, n.9 (1998) 1. 

[2] G.M.Dallavalle, Acta Phys.Pol. 28 (1997) 901 and references therein; 
B.Bloch, talk given at XXIX International Conference on High Energy 
Physics, Vancouver, July 23-29 1998. 



24 



[3] S.Jadach, O.Nicrosini et al, in Physics at LEP2, CERN 96-01 

G. Altarelli, T.Sjostrand and F.Zwirner eds., Vol.2 pg.229; 

A. Arbuzov et al., Phys.Lett.B 389 (1996) 129. 

[4] B.F.L.Ward, S.Jadach, M.Melles, S.A.Yost, pTep-ph/ 98 112451 , talk given 
by B.F.L.Ward at XXIX International Conference on High Energy 
Physics, Vancouver, July 23-29 1998. 

[5] S.Jadach, W.Placzek, B.F.L.Ward, Phys.Lett.B 353 (1995) 349; 

S.Jadach, E.Richter Was, B.F.L.Ward, Z.Was, Phys.Lett.B 353 (1995) 
362; 

S.Jadach et al, Comput.Phys.Comm. 102 (1997) 229. 

[6] S.Eidelman, F.Jegerlehner, Z.Phys.C 67 (1995) 585; 

H. Burkhardt, B.Pietrzyk, Phys.Lett.B, 356 (1995) 398. 

[7] G.Montagna, O.Nicrosini, F.Piccinini, Phys.Lett.B 385 (1996) 348. See 
also ref. 0. 

[8] S.Jadach, M.Melles, B.F.L.Ward, S.A.Yost, Phys.Lett.B 377 (1996) 168, 
Nucl.Phys.Proc.Suppl. 51C (1996) 164; 

S.Jadach et al., in Proceedings of the 28th International Conference on 
High Enrgy Physics, Z.Ajduk and A.K.Wroblewski eds. (World Sientific, 
Singapore 1997), Vol.II,pg.l072; 

B. F.L.Ward et al., Acta Phys.Pol.B 28 (1997) 925. 

[9] R.Alemany, M.Davier, A.Hocker, Eur.Phys.J '.C2 (1998) 123; 

M.Davier, A.Hocker, Phys.Lett.B 435 (1998) 427, |hep-ph/9805470| ; 



J.H.Kiihn, M.Steinhauser, MPI-PHT-98, [hep-ph/9802241 



S.Groote, J.G.Korner, K.Shilcher, N.F.Nasrallah, MZ-TH-98-02, 



hep-ph/9802374 



J.Erler, UPR-796T, |hep-ph/9803453| . 

[10] S.Jadach, M.Skrzypek, B.F.L.Ward Phys.Rev.D 55 (1997) 1206. 

[11] A.Arbuzov et al., J.Exp.Theor.Phys. 81 (1995) 638, Nucl.Phys. 
Proc.Suppl. 51C (1996) 154, Nucl.Phys. B 485 (1997) 457, Phys.Lett.B 
394 (1997) 218, Phys.Lett.B 399 (1997) 312; 
N.P.Merenkov et al., Acta Phys.Pol.B 28 (1997) 491. 

[12] S.Jadach, M.Skrzypek, B.F.L.Ward Phys.Rev.D 47 (1993) 3733. 



25 



[13] C.F.Weizsacker, Z.Phys. 88 (1934) 612; 

E. J.Williams, Phys.Rev. 45 (1934) 729. 

See also S.Frixione, M.Mangano, P.Nason, G.Ridolfi, Phys.Lett.B 319 
(1993) 339. 

[14] V.N.Baier, V.S.Fadin, V.A.Khoze, E.A.Kuraev, Phys.Rep. 78 (1981) 
293 and references therein. 

[15] F.Caravaglios, M.Moretti, Phys.Lett.B 358 (1995) 332. 

[16] F.A.Berends, P.H.Daverveldt, R.Kleiss, Nucl.Phys.B 253 (1985) 421, 
Comput.Phys.Comm. 40 (1986) 271,285,309; 

P.Aurence, G.A.Schuler et al, in Physics at LEP2, CERN 96-01 
G.Altarelli, T.Sjostrand and F.Zwirner eds., Vol.1 pg.291; 
L.Lonnblad, M.Seymour et al., ibidem, Vol.2 pg.187; 
A.Courau, talk given at VI EuroDAQNE Collaboration Meeting, 
Barcelona, November 7-9 1997. 

[17] F.Boudjema, B.Mele et al, in Physics at LEP2, CERN 96-01 
G.Altarelli, T.Sjostrand and F.Zwirner eds., Vol.1 pg.207; 
D.Bardin, R.Kleiss et al., ibidem, Vol.2 pg.3. 

[18] M.Moretti, Nucl.Phys.B 484 (1996) 3; 

F. Caravaglios, M.Moretti, Z.Phys. C 74 (1997) 291; 

G. Montagna, M.Moretti, O.Nicrosini, F.Piccinini, Eur.Phys.I.C 2 
(1998) 483; 

G.Montagna, M.Moretti, O.Nicrosini, F.Piccinini, FNT/T-98/07, 
CERN-TH/98-238, [hep-ph/9807465| , to appear in Nucl.Phys.B; 

F. Caravaglios, M.L.Mangano, M.Moretti, R.Pittau, [hep-ph/ 9807570 , 
CERN-TH/98-249, to appear in Nucl.Phys.B. 

[19] R.Barbieri, J.A.Mignaco, E.Remiddi, Nuovo Cimento 11A (1972) 824; 

G. J.H.Burgers, Phys.Lett.B 164B (1985) 167. 

[20] E.A.Kuraev, V.Fadin, Sov.J.Nucl.Phys. 41 (1985) 466; 

G.Altarelli, G.Martinelli, in Physics at LEP2, CERN 86-02 J.Ellis, 
R.Peccei eds., Vol.1 pg.47; 

O.Nicrosini, L.Trentadue, Phys.Lett.B 196 (1987) 551, Z.Phys. C 39 
(1998) 479; 

F.A.Berends, G. Burgers, W.L.van Neerven, Nucl.Phys.B 297 (1988) 



26 



429; 

S.Jadach, M.Skrzypek, Z.Phys. C 49 (1991) 577; 
M.Skrzypek, Acta Phys.Pol.B 23 (1992) 135; 

M.Cacciari, A.Deandrea, G.Montagna, O.Nicrosini Europhys.Lett. 17 
(1992) 123. 



[21] G.Passarino, |hep-ph/9810416 . 



27 



